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KCNQl channels play vital roles in cardiovascular, gastric and other systems. The conductance and 
dynamics of KCNQl could be modulated by different single transmembrane helical auxiliary proteins (such 
as KCNEl, KCNE2 and others). In this study, detail KCNQl function modulations by different regions of 
KCNEl or KCNE2 were examined using combinational methods of electrophysiology, 
immunofluorescence, solution NMR and related backbone flexibility analysis. In the presence of KCNE2 
N-terminus, decreased surface expression and consequent low activities of KCNQl were observed. The 
transmembrane domains (TMDs) of KCNEl and KCNE2 were illustrated to associate with the KCNQl 
channel in different modes: Ile64 in KCNE2-TMD interacting with Phe340 and Phe275 in KCNQl, while 
two pairs of interacting residues (Phe340-Thr58 and Ala244-Tyr65) in the KCNQl/KCNEl complex. The 
KCNEl C-terminus could modulate gating property of KCNQl, whereas KCNE2 C- terminus had only 
minimal influences on KCNQl. All of the results demonstrated different KCNQl function modulations by 
different regions of the two auxiliary proteins. 

Voltage-gated potassium KCNQl (Kv7.1 or KvLQTl) channels are widely expressed in various tissues^'^ 
such as the brain, heart, pancreas, stomach and intestine. They play crucial roles in the regulation of 
membrane potential and cellular excitability^'^. The KCNQl channel often associates with different KCNE 
auxiliary subunits to form KCNQl /KCNEs complexes, resulting in significant changes in gating properties''^. In 
cardiac cells, KCNEl associates with KCNQl to form channels having slow delayed rectifier Iks currents. The 
KCNQl /KCNE2 channel is essential for gastric acid secretion in gastric parietal cells, especially for the pH 
sensitivity and for the generation of constitutive K^ currents across the cell membrane^"^. 

Although KCNEl and KCNE2 both contain a single transmembrane helix, they may modulate KCNQl 
conductance in different ways^°. To illustrate different KCNQl modulations by KCNEl or KCNE2, structural 
comparisons of the two auxiliary proteins and electrophysiological studies of KCNQl/KCNEl or KCNQl/ 
KCNE2 complexes should be conducted. Previous electrophysiological and biochemical studies of KCNQl/ 
KCNEl interactions have shown that the transmembrane domain (TMD) and C-terminus of KCNEl might 
be crucial for association and modulation of KCNQl function^ Fluorescence resonance energy transfer 
(FRET), co-immunoprecipitation or physiology studies verified that the cytoplasmic tail of KCNEl could interact 
directly with the pore region of KCNQl or move close to the C-terminus of KCNQl ^'"^^. The availability of the 
KCNEl solution NMR structure^^'^^ enabled the docking analysis of the KCNEl-TMD to the tetrameric KCNQl 
channeP°. However, the detailed structure of KCNE2 and the structural basis of the KCNQl function modula- 
tions by KCNE2 remain unclear. Furthermore, the KCNEl solution NMR structure was determined alone in 
detergent micelles and could not explicitly explain the association mode with KCNQl in the KCNQl/KCNEl 
complex. Besides docking analysis with KCNQl channels, flexibility analysis of KCNEl or KCNE2 might also 
provide information for potential interactions with KCNQl in KCNQl /KCNEs complexes. 

In the present study, we examined the KCNEl or KCNE2 mediated regional modulation of KCNQl channel 
conductance by generating different chimeras of KCNEl and KCNE2 to clarify the roles of the N- terminal 
regions, TMDs or C-terminal tails in the association with KCNQl. Immunofluoresence data showed that the 
N-terminal region of KCNE2 (E2(N)) affected KCNQl trafficking to the cell membrane. Using Q-scanning 
mutation with double mutant cycle analysis and molecular dynamics (MD) based on the NMR structure of the 
TMDs of two KCNE subunits, models of the different association modes of KCNEl-TMD or KCNE2-TMD to 
KCNQl were developed. Comparisons of structural flexibility between the C-terminal tails of KCNEl or KCNE2 
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and their deletion mutants suggested that the C-terminal tail of 
KCNE2 plays a minor role in the function of the KCNQl channel. 

Results 

KCNE1/KCNE2 chimeras modulate the channel properties of 
KCNQl. Sequence alignment of full length KCNEl and KCNE2 
showed 32% identity and 92% similarity among the primary 
sequences (Fig. SI). The significant variation in the currents 
between KCNQl/KCNEl and KCNQ1/KCNE2 could be related to 
different modes of association between KCNEl or KCNE2 and 
KCNQl. To examine whether the subdomains of KCNEs (N 
terminus, TMD or C terminal tail) are involved in the modulation 
of KCNQl current, chimeras of KCNEl and KCNE2 were 
constructed (Fig. S2) and co-expressed with KCNQl in HEK293 
cells. 

A previous study showed that the TMD of KCNEl interacted with 
the S6 domain of KCNQP\ Two chimeras. El [E2(T)] and E2[E1(T)] 
(defined in Fig. S2a), were constructed in which only the TMD was 
exchanged to examine the possible interactions between the TMDs of 
KCNEs and KCNQl. When coexpressed with KCNQl, E2[E1(T)] 
showed slowly activated currents with attenuated amplitude, 
whereas E1[E2(T)] showed rapidly activated currents (Fig. lieft), 
indicating different KCNQl modulation by E1(T) or E2(T). 

The role of the C terminal tail of KCNEl or KCNE2 on channel 
modulation was examined with the chimeras E2[E1(C)] and 
E1[E2(C)] (defined in Fig. S2a). When coexpressed with KCNQl, 
E1[E2(C)] fully eliminated the slowly activated currents, whereas 
E2[E1(C)] showed slowly activated and high amplitude currents 
(Fig. 1 right) J indicating that the KCNEl C- terminus might modulate 
gating properties of KCNQl /KCNE channels, consistent with a pre- 
vious report^^. 

Role of the KCNE2 N-terminus in KCNQl surface expression. 

More detail electrophysiological studies of the chimera of E1[E2(N)] 
(defined in Fig. S2a), co-expressing with KCNQl, were conducted. 
The observed peak current densities were 75 ± 6.4 pA/pF (KCNQl 
alone, n = 25), 37 ± 2.8 pA/pF (Q1/E2, n = 31) and 31 ± 3.2 pA/pF 
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Figure 1 | The currents of KCNQl co-expressed with KCNEl, KCNE2 
chimeras in HEK293 cells. Representative current traces of Ql/El [E2(T)], 
Q1/E1[E2(C)],Q1/E2[E1(T)] and Q1/E2[E1(C)]. The current traces were 
elicited by voltage steps ranging from — 120 to +60 mV with 20 mV 
increments for 1.5 s (and 2 s for Q1/E2[E1(C)] due to slow activation of 
this channel) from a holding potential of —80 mV. The voltage protocols 
for these chimeras are plotted in the bottom of each trace. The Ql/ 
E2 [El (C) ] currents were evoked by a voltage protocol placed at the bottom 
in which the repolarizing voltage was — 160 mV. All the currents were 
expressed as current density in pA/pF. 



(Ql/El [E2(N)], n = 16) for the different channels at +60 mV 
(Fig. 2a, b). The experimental results showed that the Ql/El [E2(N)] 
current decreased to the background level of Q1/E2. This suggested 
that the KCNE2 N-terminus is responsible for the attenuated KCNQl 
current. 

To reveal the mechanisms of the attenuated currents of KCNQl/ 
E1[E2(N)] chimera and the suppressed currents of KCNQ1/KCNE2 
complex, surface expression of KCNQl in the presence of KCNE2 or 
chimeras were analyzed using immunofluorescence. The fluor- 
escence intensities were measured as 100 ± 4.3% (n = 194) and 
13.8 ± 1.4% (n = 198) for KCNQl and KCNQ1/KCNE2, respect- 
ively. The chimera El [E2(N)] also induced the fluorescence intensity 
reduction (15.6 ± 1.7%, n = 174, Fig. 2c, d). The relative fluorescence 
intensity ratios were I(qi/e2)/I(qi) = 0.14 and I(qi/ei[e2(n)])/I(qi) = 
0.16 (Fig. 2d). Previously, it was reported that surface expression of 
KCNQl was not affected by KCNEP^. Therefore, it is hypothesized 
that the presence of KCNE2 N-terminus might cause decreased sur- 
face expression of KCNQl in cell membrane. 

Identification of interaction sites within the KCNE-TMDs and 
KCNQl. To identif)^ the interaction sites within the KCNE-TMDs 
and KCNQl, glutamine-scanning and double mutant cycle analysis 
were conducted for all the point-mutations. Several amino acids were 
individually mutated to alanine because of low expression levels. 

In accordance with a previous report by Strutz-Seebohm et at. 
(2011), residues on the KCNEl-TMD (F57-I66) and KCNQ1-S6 
(S338-G345) were firstly scanned. The V50 values were 27 mV for 
Ql/El (n = 30), -4 mV for Q1/E1(T58Q) (n = 15), 129 mV for 
Q1(F340Q)/E1 (n = 16) and 66 mV for Q1(F340Q)/E1(T58Q) (n = 
16) (Fig. 3a). The conferred AAV50 = 32 mV (Fig. 3 A) strongly 
indicated that the residue pair of Q1(F340)-E1(T58) were involved 
in KCNQl/KCNEl interactions, which was consistent with previous 
results'". Similarly, the V50 values of 60 mV for Q1/E1(Y65Q) (n = 
13), 23 mVforQl(A344Q)/El (n = 12) and 18 mVfor Q1(A344Q)/ 
E1(Y65Q) (n = 15) conferred a AAV50 = 38 mV (Fig. 3b), indicating 
that the second pair of interaction residues was Q1(A344)-E1(Y65). 
A summary of the AAV50 values was plotted for mutant combina- 
tions in KCNQl -S6 and KCNEl-TMD (Fig. 3c). Only the two pairs 
of interaction residues [Q1(F340)-E1(T58), Q1(A344)-E1(Y65)] 
were distinctive. 

Double mutant cycle analysis was also conducted for KCNQl and 
KCNE2-TMD. The fitted V50 values of -22 mV for Q1/E2 (n = 31), 
-27 mVforQl/E2(I64Q) (n = 16), -12 mV for Q1(F340Q)/E2 (n 
= 20) and 47 mV for Q1(F340Q)/E2(I64Q) (n = 15) conferred a 
AAV50 = —64 mV (Fig. 4a), indicating an interaction pair of 
Q1(F340) and E2(I64) in the KCNQ1/KCNE2 complex. At the same 
time, the fitted V50 values of -27 mV for Q1/E2(I64Q) (n = 16), 
-82 mV for Q1(F275Q)/E2 (n = 14) and - 15 mV for Q1(F275Q)/ 
E2(I64Q) (n = 12) conferred a AAV50 = -72 mV (Fig. 4b), indi- 
cating another interaction pair of Q1(F275)-E2(I64). Therefore, the 
residue E2(I64) could interact synchronously with both Q1(F340) 
and Q1(F275). The AAV50 values were plotted for mutant combina- 
tions of KCNQl and KCNE2-TMD (Fig. 4c). Only two pairs of 
interaction residues [Q1(F275)-E2(I64), Q1(F340)-E2(I64)] within 
the KCNQl and KCNE2-TMD were mapped out using this method. 

Solution NMR structure determination of KCNE2 in detergent 
micelles. Differences between the KCNQl channel modulations by 
KCNEl or KCNE2 can be examined by analyzing their interaction 
modes, especially the spatial contacts between KCNQl and KCNEs. 
The solution NMR structure of KCNEl is available^^, and in the 
present study the solution NMR structure of KCNE2 was 
determined, providing the basis for further analysis of the different 
modes of interaction between KCNQl and the two KCNE proteins at 
the atomic structure level. 

Solution NMR structure analysis of KCNE2 in detergent micelles 
resulted in the assignment of approximately 90% of the backbone 



SCIENTIFIC REPORTS | 4 : 4973 | DOI: 1 0.1 038/srep04973 



2 



Q1/E2 



Q1/E1[E2(N)] 

60 mV 




Q1/E1[E2(N)] 



Figure 2 | Role of KCNE2 N terminus on KCNQl in the HEK293 cell, (a) Representative current traces of KCNQl, KCNQ1/KCNE2, and KCNQl with 
chimera E1[E2(N)]. (b) Bars show averaged current density of KCNQl, KCNQ1/KCNE2, and KCNQl with chimera E1[E2(N)]. Error bars 
represent SEM. < 0.001 (Student's t test), (c) Ql-146Myc (red), EGFP (a marker, green) and merged immunofluorescence images of HEK293 cells. 
Upper row: Ql; middle row: Q1/E2; lower row: Ql/El [E2(N)] . The scale bar represents 10 [im. (d) Normalized fluorescence intensities of Ql-146Myc for 
the transfections of Ql, Q1/E2 and Q1/E1[E2(N)]. Error bars represent SEM. ***P < 0.001 (Student's t test). 



resonances (Ca, CO, HN, and N)^'. Here, the structure of KCNE2 
was determined based on 582 restraints, including 240 dihedral angle 
restraints, 282 NOE restraints, and 60 hydrogen-bond restraints. The 
statistics of the quality and precision of the final 10 lowest energy 
structures were summarized in Table 82. The chemical shifts for 
KCNE2 were deposited in the Biological Magnetic Resonance Bank 
(BMRB accession number: 18817), and the coordinates were depos- 
ited in the Protein Data Bank (PDB accession code: 2M0Q). 

The overall structure of KCNE2 was shown in Fig. 5a, containing 
three a helices: an N-terminal helix (al, residue 2-42), a transmem- 
brane helix (a2, residues 45-78) and a C-terminal helix (a3, residues 
85-118). The structures of three a helices were individually conver- 
gent (Fig. 5a and Table 82), but their relative orientations were much 
dynamic due to the lack of long distance restraints, especially in the 
presence of globular shape detergent micelles^^. 

Docking and interaction models of KCNQl /KCNE-TMDs. The 

primary sequences of KCNEl and KCNE2 were listed in Supplemen- 
tary Fig. SI and comparisons of their structures were shown in 



Fig. 5b. The ot2-a3 linker and the juxta- membrane region of a3 
were more conserved than the other regions. The transmembrane 
helix of KCNE2(45-78) was 8 residues longer than that of 
KCNE1(46-71) (Fig. 5 and Supplementary Fig. SI). A kink in the 
transmembrane helix of KCNEl was present in residues around 
Gly60, whereas it was absent in the corresponding residue Ala65 in 
KCNE2, resulting in a straight helix (Fig. 5 and Supplementary Fig. 
SI). The electrostatic potential surfaces of the transmembrane helices 
of KCNEl and KCNE2 were similar (Fig. 5), both having highly 
hydrophobic amino acids (Supplementary Fig. SI). 

The solution NMR structures of the KCNEl -TMD and KCNE2- 
TMD, together with the spatial proximity restraints from the double 
mutant cycle analysis were used to construct complex models of 
KCNQl/KCNEl-TMD and KCNQ1/KCNE2-TMD using PyMol 
suite (Fig. 6a). The docked complex structure optimizations were 
processed using 2.4 ns DESMOND suite molecular dynamic simula- 
tions for KCNQl/KCNEl-TMD or KCNQl /KCNE2-TMD21 The 
optimized docking results were shown in Figure 6. In Fig. 6aieft, 
the distances between residue pairs as measured in the double 
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Figure 3 | Identification of interaction sites between KCNQl and KCNEl using double-mutant cycle analysis, (a) Left, representative current traces for 
Ql/El, Ql/El (T58Q), Ql (F340Q)/E1 and Ql (F340Q)/E1 (T58Q) as indicated. The voltage protocols are plotted at the bottom. Right, the corresponding 
GV curves from peak currents plotted for Ql/El (black), Ql/El (T58Q) (blue), Q1(F340Q)/E1 (red) and Q1(F340Q)/E1(T58Q) (green) as indicated, (b) 
Left, representative current traces for Ql/El, Q1/E1(Y65Q), Ql (A344Q)/E1 and Q1(A344Q)/E1 (Y65Q) as indicated. The voltage protocols are plotted at 
the bottom. Right, the corresponding GV curves plotted for Ql/El (black), Ql/El (Y65Q) (blue), Q1(A344Q)/E1 (red) and Q1(A344Q)/E1(Y65Q) 
(green) as indicated, (c) Identification of interaction sites for the Q1(S6) mutants (S338 to G345) with the El(TM) mutants (F57 to 166). The AAV50 
values are plotted for the mutant combinations between Ql and El as indicated. Detailed information on the V50 values is provided in Table SI. Two 
dotted lines are located on 20 and —20 mV. AAV50 values higher than 20 mV or lower than —20 mV are considered as significant changes. 



mutant cycle analysis of KCNQl /KCNE2-TMD were 3.8 A and 
2.6 A for Q1(F340)-E2(I64) and Q1(F275)-E2(I64), respectively. In 
Figure 6ai.ightJ the distances were 3.4 A and 2.9 A for Q1(F340)- 
E1(T58) and Q1(344)-E1(Y65), respectively. Interestingly, KCNEl- 
TMD and KCNE2-TMD (Fig. 6b) bound to the S6 helix of KCNQl in 
an apparently different orientation. Two restrained residues between 
KCNQl and KCNEl-TMD [Q1(F340)-E1(T58) or Q1(A344)- 
E1(Y65)] drove interactions between KCNEl-TMD and the S4-S5 
linker in a Chinese chopsticks mode (Fig. 6bj.ight)j which might con- 
strain the movement of S6 and result in the slow activation of 



KCNQl channels. On the contrary, KCNE2-TMD only provided 
one residue site (164) to interact with two sites in KCNQl (F340, 
F275), which could not restrain the movement of the S6 helix of 
KCNQl (Fig. 6bieft). 

Role of the KCNE C-terminal tails in KCNQl modulation. In the 

docking models, the KCNEl-TMD and KCNE2-TMD showed 
different association modes with KCNQl (Fig. 6b), which might 
lead to different orientations of the C-terminal tails of KCNEl and 
KCNE2. At the same time, the physiological studies between KCNQl 
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Figure 4 | Identification of interaction sites between KCNQl and KCNE2 using double-mutant cycle analysis, (a) Left, representative current traces for 
Q1/E2, Q1/E2(I64Q), Q1(F340Q)/E2 and Q1(F340Q)/E2(I64Q) as indicated. The voltage protocols are plotted at the bottom. Right, the corresponding 
GV curves from peak currents plotted for Q1/E2 (black), Q1/E2(I64Q) (blue), Q1(F340Q)/E2 (red) and Q1(F340Q)/E2(I64Q) (green) as indicated, 
(b) Left, representative current traces for Q1/E2, Q1/E2(I64Q), Q1(F275Q)/E2 and Q1(F275Q)/E2(I64Q) as indicated. The voltage protocol is plotted at 
the bottom. Right, the corresponding GV curves from peak currents plotted for Q1/E2 (black), Q1/E2(I64Q) (blue), Q1(F275A)/E2 (red) and 
Q1(F275A)/E2(I64Q) (green) as indicated, (c) Identification of interaction sites for the Q1(S6) mutants (F332, F335 and S338 to G345), the Q1(S5) 
mutants (F270 and F275) with the E2(TM) mutants (163 to V72). The AAV50 values are plotted for the mutant combinations between Ql and E2 as 
indicated. Detailed information on the V50 values is provided in Table SI. Two dotted lines are located on 20 and —20 mV. AAV50 values higher than 
20 mV or lower than —20 mV are considered as significant changes. 



and KCNE chimeras (Fig. la) indicated that the C- terminus of 
KCNEl might modulate gating properties of KCNQl /KCNE channel. 

To verify the possibilities, residue specific flexibilities of KCNE2 
proteins were analyzed using the solution NMR backbone amide^^N 
longitudinal relaxation (Tl), transverse relaxation (T2) and back- 



bone amide^H-^^N inter-nuclear dipolar coupling (NOE) measure- 
ments (Fig. 7). Compared with previously reported backbone 
relaxation data of KCNEl the C-terminus of KCNE2 was not that 
flexible. Therefore, high flexibility of the KCNEl C-terminus pro- 
vided possibilities for interaction with KCNQl. 
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Figure 5 | Solution structure of KCNE2 in detergent micelles and comparison of the structures and electrostatic potential surfaces of the KCNE 
transmembrane helical domains, (a) Solution structure of KCNE2. Left panel, Schematic illustration of the structure of KCNE2. Middle and right panels, 
The backbone superposition of the final 10 conformers with the lowest energies in the transmembrane helix (middle panel, residue 45-78), the N-terminal 
helix (right upper panel, residues 2-42), and C-terminal helix (right lower panel, 85-118) is shown. The N-terminal helix of KCNE2 is an amphiphilic 
helix. The backbone is shown in blue, and the hydrophobic and hydrophilic residues are respectively shown in red and green, (b) Comparison of the 
structures and electrostatic potential surfaces of the KCNE transmembrane helical domains KCNE2 (left panel) and KCNEl (right panel). Red 
denotes negative potential, blue denotes positive potential, and white denotes neutral potential. The side-chains of several important residues are 
indicated. 



Further examinations of different roles of the flexible KCNEl C- 
terminus and the relatively rigid KCNE2 C-terminus in KCNQl 
modulation were conducted. KCNEl with truncations in the C-ter- 
minus after residues 81, 97 and 113 was co- expressed with KCNQl to 
identify the region critical for KCNQl current regulation. The 
acquired channel currents were shown as indicated: Ql/E 1(1-81) 
(Fig. 8aieft), Ql/El(l-97) (Fig. 8a^iddie) and Q1/E1(1-113) 
(Fig. 8aj-ight)- All the constructs resulted in large slowly activated 
currents, with a significant right-shift in the GV curve (Fig. 8b). 
The increased GV right-shift with progressive truncation (Fig. 8b) 
suggested that the entire KCNEl C-terminus may be important for 
channel modulation. Channel currents were also analyzed for two 
KCNE2 truncations (E2(l-75) and E2(l-112)) coexpressed with 
KCNQl. Rapidly activated currents were observed similar as those 
elicited by wild-type KCNE2, with only minor right-shifts in the GV 
curve (Fig. 8c). Therefore, the more pronounced KCNQl kinetic 
modulations by KCNEl C-terminal truncations than KCNE2 C-ter- 
minal truncations, together with higher flexibilities of KCNEl C- 
terminus than KCNE2 C-terminus derived from solution NMR 
relaxation studies, implied that the rigid KCNE2 C-terminus might 
lead to less interaction, whereas the flexible KCNEl C-terminus 
provided a mechanism for KCNQl association and modulation. 



Discussion 

In previous reports, coassembly of KCNQl with KCNE2 can form a 
background potassium channel that is functionally distinct from the 
slow- activated high amplitude KCNQl /KCNEl channel, as 
observed in Fig. 2a^°'^^. Further studies demonstrated that association 
of KCNE2 with KCNQl/KCNEl channels could suppress the cur- 
rent amplitude^^. The observed current reduction could be attributed 
to two possible factors: (1) the fewer channel expression on the cell 
membrane due to impaired trafficking or increased internalization; 
(2) the reduced single channel conductance of each KCNQl in the 
presence of KCNE2. 

Different from the membrane targeting incapability of KCNEl, 
the KCNE2 alone can be found in cell membrane^^ When KCNE2 
associated with KCNQl, the impaired KCNE2 membrane location 
and partial targeting of KCNE2 to lipid rafts were observed^^. 
Alternatively, the decreased amount of KCNQl were detected in 
lipid rafts in the presence of KCNE2, whereas no obvious changes 
in the abundance of KCNQl in lipid rafts were observed when 
KCNQl associated with KCNEl^^. Consistent with these reports, 
KCNQl was significantly down -regulated when it co-expressed with 
KCNE2 (Fig. 2), confirming that association with KCNE2 could 
bring negative influences on the membrane targeting of KCNQl. 



SCIENTIFIC REPORTS | 4 : 4973 | DOI: 1 0.1 038/srep04973 



6 




Figure 6 | Interaction models between the transmembrane domain of KCNEs and the KCNQl channel, (a) The top panel shows the side-view of the Ql/ 
E2 (left) and Ql/El (right) complexes. The lower panel shows the top view with the distances between the inter-residues of Ql to E2 (left) and Ql to El 
(right). The sticks represent the interacting residues and the dotted lines denote distances in A. (b) Left, schematic illustration showing that the E2 has no 
site interacting with the C-terminus of S6 or the S4-S5 linker. Right, schematic illustration showing that the El interacts with S6 or the S4-S5 
linker to obstruct the channel opening via a spring-like connection. 



Furthermore, the association of KCNQl with chimera of E1[E2(N)] 
also resulted in decreased channel current and reduced expression 
level of KCNQl. These observations strongly indicated that the N- 
terminal region of KCNE2 might be responsible for the down-regu- 
lation of KCNQl membrane targeting. Detail mechanism of how the 
KCNE2 N-terminus affected the KCNQl trafficking or channel 
behaviors should be in more concerns and awaited further studies. 

Physiology studies of KCNQl associated with chimeras El [E2(T)] 
or E2[E1(T)] resulted in much different channel kinetics from the 
wild type KCNEl or KCNE2 (Fig. 1), strongly indicating that both 
KCNEl-TMD and KCNE2-TMD can associate with and modulate 
KCNQl. 

Several previous reports have demonstrated that the KCNEl- 
TMD interacts with the S6 segment or S4/S5 linker of KCNQl. 
Residues 338-340 on the S6 segment of KCNQl and residues 57- 
59 on the TMD of KCNEl were revealed as the interaction sites for 
the gating effect^°'^^'^°. The interaction between the residue pair F340- 
T58'° and residue T58 on KCNEl with residues on KCNQl S4 
domain and S5/S6 segment have been shown previously^°. In the 



present study, the results of double mutant cycle analysis showed 
that the KCNQl /KCNEl channel had two pairs of interaction sites: 
the previously reported F340-T58 pair and a new A344-Y65 pair that 
was mapped to the S6 segment of KCNQl and the TMD of KCNEl. 
As in our proposed Chinese chopstick model, the two interaction 
sites between KCNEl and KCNQl could effectively hinder the 
movement of the S6 segment, thus affecting channel activity. At 
the same time, only one residue (164, located at similar position as 
T58 of KCNEP^) was mapped on KCNE2-TMD to interact with two 
sites (F340 on S6 and F275 on S5) on KCNQl, suggesting that 
KCNE2 was positioned between the S5/S6 segment of KCNQl. 
Despite the high similarity between the primary sequences of 
KCNEl-TMD and KCNE2-TMD (Fig. SI), our docking models indi- 
cated that the straight rigid helix of KCNE2-TMD (Fig. 5) and the 
rigid KCNE2 C-terminus (Fig. 7) determined membrane crossing 
and interaction with KCNQl without physical interaction with cyto- 
plasmic regions of KCNQl. Consequently, KCNE2 has a lesser effect 
on KCNQl channel gating. On the contrary, the kinked helix in 
KCNEl-TMD (Fig. 5) and the flexible KCNEl C-terminus'" 
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Figure 7 | Backbone relaxation measurement of the KCNE2 protein in detergent micelles. Backbone ^^N Tl longitudinal relaxation and T2 
transverse relaxation, stable ^H-^^N NOE effect were plotted for each residue of the protein. The transmembrane helical domain is delimited by two dotted 
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result in a tighter interaction with KCNQl (two interaction pairs) 
and direct contact with the cytosolic region of KCNQl, thus slowing 
channel activation. 

Moreover, the residues 67-78 in KCNEl C-terminus were 
mapped to be important for KCNQl gating modulation^^. Deletion 
truncations after residue 70 were observed not to alter the activation 
profile of Iksj but significantly right-shifted the voltage dependence 
of the channel activation by 60 mV^^. Physiology studies of KCNQl 
associated with chimera E2[E1(C)] induced a slow activated KCNQl 
current, while data of KCNQl associated with chimera E1[E2(C)] 
retained the fast activation conductance of KCNQl. These observa- 
tions indicating that the C-terminus of KCNEl can interact directly 
with the KCNQl channel. 

Residue specific relaxation analysis of KCNEl and KCNE2 in 
detergent micells using solution NMR demonstrated obvious differ- 
ent dynamics of C- terminal tails of the two channel auxiliary pro- 
teins: high rigidity in KCNE2 C-terminus versus flexibility in KCNEl 
C-terminus. Biophysically, the flexible KCNEl C-terminus provided 
possibilities for its interaction with some segments of KCNQl and 
consequent KCNQl function modulations. To verify this hypothesis, 
physiological studies of KCNQl associated with KCNEs with differ- 
ent C-terminal truncations were conducted. The activation kinetics 
of the KCNQl /KCNEl channel remained in a slow activated current 
with all three truncations El(l-81), El(l-97) and El(l-113) 
(Fig. 8a), whereas the KCNEl truncations El(l-81) and El(l-97) 



shifted the G-V curves of the activation currents to more positive 
potentials by approximately 60 mV, similar as the previously 
reported KCNEl deletion mutants after residue 70 By contrast, 
the G-V curves of the KCNE2 C-terminus truncations E2(l-75) 
and E2(l-112) showed only a slight right shift. 

In a summary, chimeras were generated and Q-scanning methods 
were used to understand the mechanism of KCNQl current regu- 
lation by KCNEs. The results of physiological analysis suggested 
that the KCNE2 N-terminus was responsible for the decreased 
amplitude of the KCNQ1/KCNE2 channel, and the TMD and C 
terminus of KCNEl were involved in the regulation of KCNQl/ 
KCNEl channel kinetics. Comparisons of solution NMR structures 
of KCNEl and KCNE2 showed a straight helix in the KCNE2-TMD, 
whereas the KCNEl -TMD had a kinked helix. Modeling and dock- 
ing analysis suggested different association modes of KCNQl/ 
KCNE2 versus KCNQl/KCNEl, leading to the differential modu- 
lation of KCNQl channel gating. Solution NMR backbone relaxa- 
tion analysis of KCNE2 demonstrated a more rigid C-terminal tail, 
which was different from previous reported flexible KCNEl C-ter- 
minus. Truncation and physiological studies verified the essential 
roles of the flexible KCNEl C-terminus in KCNQl modulations. 
Although different regions of KCNEl or KCNE2 were mapped for 
KCNQl modulations, the detail mechanisms on how these seg- 
ments affected KCNQl conductance and kinetics require further 
studies. 
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Figure 8 | Currents and G-V curves for Ql /El -truncations, Q1/E2, and 
Ql/E2-truncations. (a) The currents of Q1/E1(1-81) (left) and Q1/E1(1- 
97) (middle) were evoked by voltage steps ranging from —80 to 140 mV 
with 20-mV increments from a holding potential of — 80 mV, respectively, 
and the current of Q1/E1( 1-113) (right) was evoked by voltage steps 
ranging from — 120 to 100 mV with 20-mV increments from a holding 
potential of — 80 mV. The voltage protocols are plotted at the bottom, (b) 
The G-V curves were plotted for Ql/El, Q1/E1(1-81), Ql/El(l-97) and 
Q1/E1(1-113) as indicated, (c) The currents of Q1/E2, Ql/E2(l-75) and 
Ql/E2(l-112) were evoked by voltage steps ranging from —80 to 60 mV 
with 20-mV increments from a holding potential of —80 mV. The G-V 
curves were plotted for the three complexes at the bottom. 

Methods 

Chimeras and mutations. Constructs of the wild-type human KCNQl subunit or 
KCNQl-146Myc were prepared as described previously^^. The coding sequences of 
human kcnel and kcne2 (GenBank accession NOs NG_009091 and NG_008804) 
were sub-cloned into pcDNA3.1( + ) (Invitrogen, CA) using the appropriate 



restriction enzymes. All chimeras were constructed using overlap -extension PCR 
with Pfu polymerase (Fermentas). Mutations were introduced using the Site-directed 
Gene Mutagenesis Kit (Beyotime). All chimeras and point mutations were verified 
through direct DNA sequencing. 

Cell culture and transient transfections in HEK293 cells. HEK293 cells were 
cultured in DMEM supplemented with 10% fetal bovine serum (FBS), 100 U/ml 
penicillin and streptomycin in 37°C and 5% CO2 incubators. One day before 
transfection, cells were transferred to 24-well plates, and when reaching 90% 
confluence, transfections were performed with Lipofectamine 2000 (Invitrogen, CA). 
At 4-6 h after transfection, cells were transferred to poly-D-lysine (Sigma, MO) 
coated slides for patch-clamp recording or imaging. For all co-transfection 
experiments, the ratio of KCNQl coding plasmids to KCNEl (or KCNE2) coding 
plasmids was 1 : 1.5. 

Immunofluorescence imaging in HEK293 cells. Approximately 24 h after 
transfection, cells were fixed with 2% paraformaldehyde in Phosphate Buffered Saline 
(PBS) for 10 min, washed with PBS (5 min X 3), and blocked with 5% FBS for 1 h. 
Then, cells were incubated with mouse monoclonal anti-human Myc antibody 
(1 : 300) (Abeam) overnight at 4°C, washed with 2.5% FBS in PBS (5 min X 6), and 
incubated with Rhodamine-conjugated goat anti-mouse IgG (H + L) (1 : 300) 
(Proteintech Group) for 2 h at room temperature (25°C). Cells were then washed 
(5 min X 6) and soaked in PBS, then were visualized using a confocal laser scanning 
microscopy (Olympus 1X71) and a lOOX oil immersion lens (NA1.30). Parameter 
selection, sample scanning and image acquisition were all controlled by the Andor IQ 
2.0 software (Andor Technology pic, Springs^ale Business Park, United Kingdom). 16- 
bit digital images were obtained with an EM CCD camera (DU897K, ANDOR iXon, 
United Kingdom). The quantification of the average fluorescence intensity of 
KCNQl -Myc was performed in ImageJ 1.45 (Wayne Rasband, National Institutes of 
Health). Following background subtraction, a 4-pixel- wide region along plasma 
membrane in the middle layer was drawn for the quantification. Adjacent membranes 
from two transfected cells were rejected from the analysis. Fields for imaging and 
analysis were selected using EGFP fluorescence to identify transfected cells. All 
experiments were performed at least three times. For each group, approximately 60 
cells were examined and normalized to the intensity of the control group. 
Representative images were deconvolved with AutoQuant X2 software (Media 
Cybernetics Inc., MD) and arranged by Adobe Photoshop CS5. 

Solutions. For the whole-cell recordings, the pipette solution contained the following 
(in mM): 5 Ethylene glycol-bis(2-aminoethylether)-N,N,N',N'- tetraacetic acid 
(EGTA), 140 KCl, 10 NaCl, 1 MgCb, 10 HEPES (pH 7.0) titrated with KOH. Bathing 
solution contained the following (in mM): 150 NaCl, 4 KCl, 2 CaClz, 1 MgClz, 10 
HEPES (pH 7.0) titrated with NaOH. All chemicals were purchased from Sigma. 

Patch clamp recording. Transfected HEK293 cells were bathed in a solution 
containing the following (in mM): 150 NaCl, 4 KCl, 2 CaCl2, 1 MgCla, 10 HEPES one 
day after transfection. All the experiments were carried out in whole-cell patches. 
Patch pipettes were pulled from borosilicate glass capillaries with a resistance of 2-3 
megohms and filled with the pipette solution. Experiments were performed using an 
Axon 200B patch-clamp amplifier with its software (Axon Instruments, CA, USA). 
Currents were typically digitized at 10 kHz and filtered at 5 kHz. All experiments 
were performed at room temperature (22-24°C). 

Data analysis. Recording data were analyzed with Clampfit (Axon Instruments) and 
Sigmaplot (SPSS) software. Unless stated otherwise, the data are presented as mean ± 
SE and significance was tested by Student's T test. Differences between the mean 
values were considered significant at a probability <0.05. The G-V curves for 
activation were fitted by the Boltzmann equation: G/Gmax = (1 + exp((V — V50)/ 
k))~\ where V50 is the voltage at which the conductance (G) is half the maximum 
conductance (Gmax) and k is a factor affecting the steepness of the activation. 

For double mutant cycle analysis^^"^^ the principle is that pairwise coupling free 
energy changes associated with mutations in pairs of residues located in different 
subunits can be calculated using a thermodynamic square comprised of the wild-type 
(WT) complex (aP), the two single mutants (a*P and ap*) and the corresponding 
double mutant (a*P*). Here, a, P and * denote the KCNQl, KCNE and mutants, 
respectively. The thermodynamic square was described as AAG= AG^p* + AGo^^p — 
AGa*p*. In this cycle, if the total free energy change A AG is significant (> 1 kcal/mol), 
the pairwise residues on a and P subunits can be considered as coupled^^. Since the 
free energy of a channel is directly proportional to the potential at which the open 
probability is 50% (V50); therefore, the V50 was measured and used for double mutant 
cycle analysis as follows: AAV50 = AVso-ocp* + AVso-oc^p ~ AV5o-a*p* ~ AVso-ap- In 
this method, significant changes of AAV50 ^ 20 mV (or 1 kcal/mol) can be used as 
the criteria to examine whether the two sites in KCNQl or KCNEs are involved in the 
direct interaction. 

Structure determination of KCNE2 using solution NMR. The human KCNE2 
protein was over-expressed using an E. coli bacterial expression system, and purified 
in LMPG micelles. A series of TROSY-based multi-dimensional NMR experiments 
were recorded on a 600 MHz Bruker spectrometer. The 2D^H,^^N-HSQC, 3D 
HNCO, HNCOCA, HNCA, HNCACB, and HN(CO)CACB spectra were collected 
using a'H/i'C/^'N labeled KCNE2 sample in LMPG micelles at pH 6.0 and 40°C. The 
3D HBHA(CO)NH and 3D^='N-edited NOESY (mixing time: 100 ms) spectra were 
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collected using a^^C/^^N labeled KCNE2 sample. Details of the experimental 
parameters were described previously^^. All of the NMR spectra were processed using 
NMRPipe^*^ and analyzed using NMRView^^. Backbone resonance assignments were 
achieved^^ and deposited in the BioMagResBank (www.bmrb.wisc.edu) with 
accession number 18817. Backbone structural restraints of KCNE2 in detergent 
micelles were collected in the form of the backbone resonance chemical shift and 
short-range Nuclear Overhauser Effects (NOE). After chemical shift index and 
protein backbone dihedral angle analysis using TALOS+^*^, XPLOR-NIH^'^ was 
applied to calculate the three-dimensional structure of KCNE2 in detergent micelles. 
In the 10 lowest energy structures of the 100 calculated structures, the backbone 
structures of the KCNE2 transmembrane domain [residues 46-78] were well 
converged to an RMSD value of 0.909 A, although the backbone structures of the full 
length KCNE2 were not converged because of lack of long distance restraints. The 
statistics of the determined KCNE2 structure are provided in Table S 1 . Coordinates of 
the KCNE2 ensemble were deposited in the Protein Data Bank (PDB accession code: 
2M0Q). 

Molecular modeling and MD simulations. The KCNQl homology model was 
generated using MODELLER*". Transmembrane domains (TMD) of KCNEl or 
KCNE2 obtained from NMR studies were optimized in a POPC lipid bilayer 
membrane using DESMOND'l Docking of the optimized KCNEl -TMD and 
KCNE2-TMD with the adjacent- KCNQl subunits homology model was performed 
using Rosetta*\ The KCNQl/KCNEl-TMD and KCNQ1/KCNE2-TMD were 
inserted into two different explicit POPC lipid bilayer membranes (300 K) and then 
the whole system was neutralized by addition of Cl~ ions. Molecular dynamics 
simulations were performed with DESMOND based on a standard Dynamics 
Cascade pipeline for the simulation of macromolecular systems. The detailed 
structural conformations of KCNQl/KCNEl-TMD and KCNQ1/KCNE2-TMD 
were prepared using the latest PyMOL Suite (Version 1.5.0.4 Schrodinger, LLC). 
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